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In the framework of the diagrammatic method with self-consistent 
field, the maximum on the temperature dependence of the suscep- 
tibility of a weakly doped narrow-band Hubbard magnet below the 
Curie temperature Tq is predicted. By numerical calculations, it 
is proved that it appears at small hole concentrations n^. In this 
case, the temperature dependence of the magnetization M(T) has 
a typical Fermi-like shape with the point of inflection at a temper- 
ature Tinf . The approximate solution of the system of equations 
for the mean spin and the chemical potential gives the Schottky 
susceptibility typical of a two-level system with a gap of order of 
n^(l — n^)iy, where W is the bandwidth. This behavior reflects 
the existence of two subbands with up- and down-spins. It may 
be observed experimentally in the surface layers of oxide metal 
nanoparticles with narrow bands and a weak oxygen nonstoichiom- 
etry. 



At present, the important role of nanoparticle surfaces 
in the magnetic ordering of the electron subsystem with 
a narrow band is well established. In particular, a model 
based on the ideal inner core and an outer shell of 
nanoparticles with oxygen nonstoichiometry, vacancies, 
and stresses was suggested for ferromagnetic metallic 
compound La2/3Cai/3Mn03 in work [1]. A weak sur- 
face ferromagnetism is a universal peculiarity of many 
oxide nanoparticles including nonmagnetic AI2O3, ZnO, 
and Ce02 [2]. A complicated character of the interac- 
tion between surface and volume magnetic ions in some 
cases causes the exchange bias in nanostructures [3]. 

The increase of the susceptibility x{T) at lowering the 
temperature T, a sharp drop of x(T) with increase in 
a magnetic field, and the maximum in the curve x(T) 
at a definite doping parameter x and a temperature 
T = Tiiiax were observed earlier in HTSC compounds 
RBa2Cu30fi+rf (where R = Y, Nd, Sm) in the dielectric 
phase [4 7]. In work [8], the experimental results were 
interpreted, by considering the competition of the in- 
plane antiferromagnetic (AFM) and the weak ferromag- 
netic (FM) interplanar exchange, whose origin is related 
to the interaction between the quadrupole center and 



virtual magnons of the AFM matrix. The maximum of 
x(T) was established [9] to occur in strongly correlated 
electron systems as a result of correlations between lo- 
calized spins in the AFM phase. With increase in x, the 
maximum shifts to the area of lower temperatures, which 
points out the strengthening of a destructive influence of 
the hole dynamics on the AFM order. 

To explain the phenomena observed in RBa2Cu306+d 
compounds, we will suppose the existence of a heteroge- 
neous magnetic structure in grains of the ceramics under 
study. This magnetic heterogeneity is seemingly related 
to the structural inhomogeneity caused by the influence 
of a grain boundary. As a result, the oxygen index can 
vary at moving from the grain center to the boundary. 
Apparently, this gives the excess of the hole content on 
the surface. Therefore, in spite of the fact that the rest 
volume has AFM order, it is destroyed near the grain 
surface with forming a weak FM or paramagnetic state 
(PM). The realization of a FM or PM state depends on 
the degree of doping and the relation between the kinetic 
energy of electrons and their indirect exchange interac- 
tion [10-12]. The comparison of calculated and experi- 
mental values of x(^) will allow us to make a conclusion 
about the surface ferromagnetism in RBa2Cu306+d- 

For the description of the low-temperature thermody- 
namics of electrons with a narrow band, the theory of 
effective field [11-13] was proposed. By means of the 
scattering matrix formalism in the first order of pertur- 
bation theory with the formal parameter t (hopping inte- 
gral), one can find correctly the kinematic contribution 
to the magnetization and the susceptibility of a weakly 
doped Hubbard magnet. The system of equations de- 
scribing the mean spin (S*^) and the chemical potential 
/U as functions of T, the bandwidth W, and the electron 
concentration n is as follows: 

(i^-O) = {F^O), - 1 ^ f(E^^) + /(£,). (1) 

Here, the mean combined population {F'^°) = 1 — 2. _)_ 
(t{S^), cr=±l, /i is a magnetic field, = —/j. — \(jh — 
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Fig. 1. Temperature dependences of the mean electron spin at 
K=0 and n = 0.8, 0.95 and 0.99 (solid curves 1-3 correspond to 
Eq. (3)). The dotted curves were obtained from Eqs. (1) by 
numerical calculations 

Kt{0) {f — a{S^)}, K is the exchange parameter, / is the 
Fermi distribution function, and 

Here, t{q) is the Fourier transform of the hopping inte- 
gral. The "dressed" population {F^°)i can be expressed 
as 



(2) 



where E„ = -e^ + (5/i^^ and Sfi^ == jfJ2 Pt{q)f{Eq„)- 

q 

In the limit of low temperatures T <C Tc, where Tc is 
the Curie temperature, the system of equations (1) was 
solved analytically [12]. In particular, we obtain 



n 
2 



(3) 



where e = Wn{l — n)/2 + h — Wk{S^.) It follows from 
Eq. (3) that a deviation (S^) from the saturated mag- 
netization obeys the Fermi distribution with the param- 
eter s characterizing a splitting of the bands with up- 
and down-spins. This splitting is formed by the applied 
magnetic and offcctivo kinematic fields competing with 
a self-consistent field of the exchange origin. 

At h=0 and k = 0, the parameter i depends on the 
factor n(l — n), whose value may be small at n -C 1. 
In this case, the thermal fluctuations can have a pro- 
nounced effect on the Nagaoka's state. In particular, the 
temperature dependence of the magnetization changes 



from the typical Fermi-like behavior with the inflection 
point at Tmax to the Brillouin-type behavior without it. 
From (3), it is easy to find 



X{T) 



dh 



1 

4T 



cosh 



Wn{l -n) + 2h 
4T 



(4) 



where k = 0. Susceptibility (4) has a Schottky-like max- 
imum at T = Tmax determined from the equation 



JW = 0.5n(l - n) tanh 1 ^^!^ 



n) 



(5) 



Unfortunately, the approximate solution (3) contains 
the inflection point always because it works at low tem- 
peratures. The exact numerical solution of Eqs. (1) at 
various concentrations n shows a change of the shape 
of (S^) as a function of T at the deflnite concentration 
n (Fig. 1). To understand the reason for such a mod- 
ification, we will obtain analytical expressions for the 
longitudinal molar susceptibility X||(T) from (1). The 
differentiation of (1) with respect to T gives 



XII (T) = 



kB 



(6) 



(7) 



where /xb is the Bohr magneton, g is the ^-factor along 
the z axis, Np^ is the Avogadro constant, and fce is Boltz- 
mann's constant. We also have 

= a/3(F,_,(l - + F„{1 - (3^^)+ 

+A_<, -/(£_,) (1 - /(£_<,))). 

Here, 

^. = iEW«.) (!-/«.)), 
V'. = iEt'(?)/«) (!-/«)) 

Q 



(8) 



and 



(9) 



2 ' 



F„ 



2 ' 



(10) 
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Fig. 2. Concentration dependences of Tc (curve 1) and Tmax 
(curve 2) 

where P = l/T. 

In the PM phase {S^) = 0, and functions (8)-(10) do 

not depend on the index a: (pa = (f, ijja = V-": = 

= e, E^ = E, F„ = F, = F^- Therefore, Eq. 

(6) is simpHfied to 
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Fig. 3. Temperature dependences of the molar susceptibility X|| (T) 
at K=0 and concentrations n = 0.8 (o), 0.95 (6), and 0.99 (c) 



{F + Fn){T -^)+TA- Tf{e){l - fje)) 

T^-i}{F + F^{) + T^ • ^ ' 

The equaUty of the denominator in (11) to zero yields 
the following expression for the Curie temperature: 



rp2 



2 + e" 



N 



Here, Eq 



(12) 



(1 - f ) t((z) and (5^0 = ^ E i{<l)!{Eq) 



Setting (S"^) = in (1), we can obtain the chemical 
potential /x as a solution of the equation 



1 



1 + 2e T 



-^f{E,) + 2f{-i,). 



Expanding the denominator of X||(^) in ^ series in 
T — Tc, one can find the Curie-Weiss law 



XII (2^) = 



c 



T-Tc' 



(13) 



where the Curie constant has the form 



c = 



kB 



{F + Fn) (Tc -^)+ TcA - Tc/(-M)(l - /(-/x)) 



Q(Tc) 



(14) 



In Fig. 1, the comparison of the values of mean spin 
calculated from the approximate equation (3) (hnes) and 
the exact equation (1) (points) was made. It follows from 
Fig. 1 that the character of the temperature dependence 
of {S^} is changed with increase in the electron concen- 
tration. 

Using Eq. (12), one can find the concentration de- 
pendence of Tc presented in Fig. 2. Here, the function 
Tmax{n) obtained by solving Eq. (5) is displayed. From 
Fig. 2, it is seen that r,„ax > Tc at n<0.88. Therefore, 
the Schottky-like anomaly in X||(T) disappears at these 
concentrations. This claim is supported by the results 
of numerical calculations of (11). 

Figure 3 shows that, at n=0.8, the maximum on curve 
X|| (T) is absent. Moreover, the absolute values of sus- 
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Fig. 4. Concentration dependence of the Curie constant C (curve 

1). Straight line 2 corresponds to 3^^^f~^ for localized electron 
spins 

ceptibility (11) exceed essentially the experimentally ob- 
served ones for SmBa2Cu306+(j [7]. This indicates that 
the volume, in which a weak ferromagnetism is observed, 
is small. 

In Fig. 4, the concentration dependence of the Curie 
constant is presented. One can see that, at n=l, we 

have the evident result C = ;|^-^f^— ^ corresponding to 
localized electron spins, and C as a function of n has non- 
monotone character. The mentioned dependence can be 
used to find the electron concentration knowing C from 

experiments. 

Hence, the temperature dependence of the longitudi- 
nal susceptibility of a Hubbard magnet at low tempera- 
tures in the limit of a weak doping shows the Schottky- 
like anomaly. We suppose that the condition for the 
appearance of this anomaly is related to a small value of 
the mean thermal energy of electron excitation from the 
down- to spin- up subband. If this energy is comparable 
with the mean electron kinetic energy of FM ordering, 
the indicated anomaly disappears. 
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HHSbKOTEMnEPATyPHHH OEPOMArHETHSM 
y CJIABKO JlErOBAHOMY XABBAPflOBCbKOMY 
MAPHETHKY 

E.e. 3y6o6 
P e 3 lo M e 

y Meacax Teopii' e4>eKTHBHoro caMoy3rofl}KeHoro nojia fljia 
cjiafiKo jieroBaHoro xa66apflOBCbKoro MarneTHKa nepeflGaieno 
noHBy HHSbKOTCMnepaTypHOro MaKCHMyiviy Tiiny lUoTTKi Ha 
KpHBifi cnpiiHHHTjiHBOCTi. rioKasaHO, Lu,o A^ny anoMajiiio 
syMOBjieHO KOHKypcHuieio KineMaTHMKHx i TenjiOBHx BsaeMO^iii 
KOjieKTHBisoBaHHx ejieKTpoHiB. On;iHKa BejmHHHH Mojiapnoi' 

CnpHHHaTJIHBOCTi fl03B0Jiae 3P06hTH BHCHOBOK npO MO>KJIHBy 

peajii3aniio 3a3HaHeHoro aBnma b noBepxHeBHx mapax 

MeTajIOOKCHflHHX HaHOHaCTHHOK Y3 BySbKHMH 30HaMH. 
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